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I. INTRODUCTION 

In addition to direct searches for new physics (NP) per¬ 
formed at very high energies at the LHC, low-energy pre¬ 
cision flavour observables provide a complementary win¬ 
dow to physics beyond the Standard Model (SM). Here, 
the anomalous magnetic moment of the muon (g — 2) M is 
a prime example as it is very sensitive to physics beyond 
the SM entering via quantum corrections. Also, tauonic 
B meson decays and r -A pvv are excellent probes of NP: 
they test lepton flavour universality - satisfied in the SM 
- and are sensitive to new degrees of freedom that couple 
proportional to the mass of the involved particles (e.g. 
Higgs bosons [1]) because of the heavy r lepton involved. 
The observation of flavour violating decays of the SM 
Higgs, most importantly h —> pr, would prove the exis¬ 
tence of physics beyond the SM. 

Let us briefly review the current experimental and the¬ 
oretical situation in these processes. The world average 
of the measurement of = (g — 2)^/2 is completely 
dominated by the Brookhaven experiment E821 [2] and 
is given by [3] 

a® xp = (116 592 091 ± 54 ± 33) x 1CT 11 , (1) 
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where the first error is statistical and the second system¬ 
atic. The current SM prediction is 

a® M = (116 591855 ± 59) x 10 -11 , (2) 

where almost the whole uncertainty is due to hadronic 
effects. 1 This amounts to a discrepancy between the SM 
and experimental value of a® xp —a® M = (236±87) x 10 —11 , 
i.e. a 2.7er deviation. It is not yet clear if this discrep¬ 
ancy is due to NP or rather underestimated hadronic 
uncertainties; there are ongoing efforts to reduce the 
model dependence in the hadronic light-by-light estimate 
based on dispersion relations [13-15] or lattice QCD [16- 
19]. Possible NP explanations besides supersymmetry 
(see for example Ref. [20] for a review) include very 
light Z' bosons [21-28], leptoquarks [29, 30], additional 
fermions [31] but also new scalar contributions in two- 
Higgs-doublet models (2HDM) [32, 33], also within the 
lepton-specific 2HDM [34-36]. 

Concerning tauonic B decays, the BaBar collabora¬ 
tion performed an analysis of the semileptonic B decays 


1 The SM prediction includes the QED corrections [4], an electro- 
weak contribution [5—7], leading order, next-to-leading order [8, 
9] and next-to-next-to leading order [10] hadronic vacuum polar¬ 
isation contributions and the model-dependent leading order [11] 
and next-to leading order [12] hadronic light-by-light contribu¬ 
tions. 
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B —> Dtu and B —> D*tv using the full available data 
set [37]. Recently, these decays have also been reana¬ 
lyzed by Belle [38, 39], and LHCb measured the mode 
B -A D* tv [40]. These experiments find for the ratios 


(17.39 ±0.04)% [66], alleviating the deviation to approx¬ 
imately 2er: 

a h fag = (o.59 ± 0.32)% . (9) 




BR (B -A TV) 

BR(R -A D(*)tv) ’ 


where i = e or fi, the following values: 


^(D)babar = 0.440 ± 0.058 ± 0.042 , 
^(D)belle = 0.375 ±0.064 ±0.026, 
■K(£>*)babar = 0.332 ±0.024 ±0.018, 
K(D*) BELLE — 0.293 ±0.038 ±0.015, 
^P^LHCb = 0.336 ± 0.027 ±0.030. 


( 3 ) 


( 4 ) 


Here, the first error is statistical and the second one 
is systematic. Combining these measurements one 
finds [41] 

K(D) exp = 0.388 ±0.047, 

H{D*) exp = 0.321 ± 0.021. 

Comparing the experimental values to the SM predic¬ 
tion 2 


K(D) SM = 0.297 ±0.017, 
K(D*) SM = 0.252 ±0.003, 


we see that there is a discrepancy of 1.8<r for 1Z(D) and 
3.3cr for 1Z(D *); combining them gives a 3.8<r deviation 
from the SM (compared to 3.4u taking into account the 
BaBar results only [49]). Models solving the TZ(D) puzzle 
have been discussed extensively in the literature [50-64], 
including the possibility of charged Higgs particles [51, 
53, 54], 

For r -A fivv, the dominant uncertainty in the SM 
prediction for the branching ratio comes from the r 
lifetime r T . Using the PDG [3] values for r lifetime, 
t t = (290.3 ± 0.5) x 10 -15 s, and branching ratios 

B p = BR(r -A A^)ex P = (17.41 ± 0.04)% , 

B e = BR(r -A eVv) exp = (17.83 ± 0.04)% , 

we can determine the deviations = Bt/Bf M — 1 [65] 
from the SM prediction to be 

A™ 0 = (0.69 ± 0.29)%, A e = (0.28 ±0.28)%. (8) 

There is a deviation of about 2.4tr in the muon data, 
whereas the electron channel is compatible with the SM 
prediction. HFAG finds essentially the same value for 
the tau lifetime, but a slightly lower BR(r —> /im/) exp = 


2 For these theory predictions we again used the updated results 
of [37], which rely on the calculations of Refs. [42, 43] based on 
the previous results of Refs. [44-48]. 


Again, charged Higgses [65, 67, 68] but also neutral Z' 
bosons [69] affect this decay. 

Moving from lepton non-universality to outright lepton 
flavour violation, we are drawn to the recent CMS excess 
of 2.4cr in h -A /rr [70]: 

BR (h -A fiT ) = (0.84±g;||) % . (10) 

Possible explanations naturally require an extended 
Higgs sector [33, 71 76]. 

As we see from the previous discussion, all discrepan¬ 
cies outlined above can be solved by additional scalar 
bosons. In the simplest extension of the SM with addi¬ 
tional charged Higgses, a 2HDM, one introduces a sec¬ 
ond Higgs doublet and obtains four additional physical 
Higgs particles (in the case of a CP conserving Higgs 
potential): the neutral CP-even Higgs H , a neutral CP- 
odd Higgs A and the two charged Higgses H ± . 2HDMs 
have been studied for many years with focus on the 
type-II models [77-79]. However, there are also other 
models without flavour-changing neutral currents at tree- 
level, i.e. type-I, type-X (lepton-specific) and the flipped 
2HDM (see Ref. [80] for an review). More general mod¬ 
els with flavour-changing neutral Higgs couplings at tree- 
level are named type-III models. Here the focus has been 
on minimal flavour violation [81-83], alignment [84, 85] 
or natural flavour conservation [82, 86] but also generic 
2HDMs of type III have been studied [54]. In these mod¬ 
els type II is recovered in the absence of flavour-changing 
neutral Higgs couplings. 3 While the type-II model can¬ 
not explain 7 Z(D) and 1Z(D*) simultaneously [37], this 
can be achieved by supplementing the model with ad¬ 
ditional non-holomorphic couplings [51]. However, this 
model (as the normal type-II model) is under pressure 
from b —> sy data [89] and LHC searches for A —> tt [90]. 
Furthermore, no sizable effect in the anomalous magnetic 
moment of the muon can be generated [54] and also ex¬ 
plaining h —> jjT is challenging. 

Therefore, we choose to consider in this article the 
2HDM of type X (lepton-specific). In this model the 
heavy Higgs couplings to quarks are suppressed com¬ 
pared to the type-II model and the bounds from b —> sy 
and LHC searches are much weaker, leaving more space 
for effects in (g— 2) M and tauonic B decays. Furthermore, 
the sign of the coupling of heavy Higgses to the r lepton 


3 The decoupling limit of the MSSM at tree-level is the 2HDM of 
type II. However, non-decoupling 1-loop corrections involving the 
Higgsino mass parameter ^ or non-holomorphic A-terms generate 
“wrong” Higgs Yukawa couplings giving rise to flavour chaning 
neutral Higgs couplings (see for example Ref. [87] for a complete 
one-loop analysis and Ref. [88] for the 2-loop SQCD corrections 
in the MSSM with generic SUSY breaking terms). 
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can be reversed, allowing for constructive interference in 
r —> Also, large effects in h —> /rr compared to the 

type-II-like model are possible as the Barr-Zee contribu¬ 
tions to r —> p -7 involving quarks are suppressed. 

The article is structured as follows: in the next section 
we outline our model, i.e. the 2HDM of type X with addi¬ 
tional non-standard Yukawa couplings. Sec. 3 discusses 
the relevant observables and collects the necessary for¬ 
mulae. Sec. 4 contains the numerical analysis. Finally 
we conclude in Sec. 5. 


II. 2HDM-X 

We will study a lepton-specific 2HDM (2HDM-X), de¬ 
fined by the Yukawa couplings in the Lagrangian 

C Y = -Q l Y u $ 2 u r - Q L Y d <$> 2 d R - L L Y e ^e R + h.c., 

( 11 ) 

with additional couplings that break the type-X structure 

A Cy = —Ql£ u ®i u r ~ QA d ®id R — L R ^^ 2 e R + h.c. 

( 12 ) 


where the couplings are given by 


Y'hLR ^ 
1 QiQj ~ 

-M 

m qi 

V 

S^ cos a — 

e-j sin , 

(16) 

-pHLR ^ 
1 QiQj ~ 

-M 

m qi 

V 

- Sij sin a + 

e ij cos a) , 

(17) 

taALR 

1 didj — 

1 

“vf 

d 

ij 5 



(18) 

-pALR ^ 

4 UiUj 

i— e u . 

V 2 13 

5 



(19) 

1 H+LR ^ 

Uid-i 

v w 4 i 

? 



( 20 ) 


RL ^ t; 

1 mdj — v j'j > 

YhLR _ sin a tan £ (m ti 

T ~ 72 \— Sfi ~ e N 

J.HLR _ cos a tan/3 (m ti t 

- 71 Kv fi ~ efi 

yAlr tan/3 (m e 

Y^ R ^ tan/3 ~ £/,:) > 


fi ’ 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


For arbitrary matrices Y u ' d,e and this simply 

parametrizes the 2HDM with the most general Yukawa 
interactions (type III). We will however assume the cou¬ 
pling structure to be close to type X, i.e. the £ matrices 
to be small perturbations. After electroweak symmetry 
breaking the following field redefinitions are necessary in 
order to render the fermion mass matrices diagonal 

dL,R -A D' l r <1l,r , 

ul,r-> U' l r ul, R , (13) 

#l,r L\,r£l,r ■ 

We define the (non-diagonal) coupling matrices 

e u = UlCU R , e d = D\i d D R , e e = L[^L R , (14) 


and express the Yukawa couplings in terms of the physical 
masses and the couplings e . 4 Note that we are not con¬ 
cerned with the issue of neutrino masses in this article, 
and hence do not introduce e.g. right-handed neutrinos. 
Therefore, the Higgs interactions with fermions can be 
written as 

C D v l r^ ] LR p R £ j H+ 

+ Ui (r ”X RLp L + T uX LRp r) d 3 H+ 

+ £ £ (fi T M J LRpR h H °) +h - C -' 


4 Note that since we eliminate the Yukawa couplings and not the 
couplings £ from the Lagrangian, the 2HDM-X is recovered in 
the limit e —> 0. 


in the limit of large tan f3 of interest in this article. 
V = U^D r denotes the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix and v ~ 174 GeV the vacuum ex¬ 
pectation value. Note that for e | 3 > m T /v the sign of the 
couplings of A , H and H + to taus is reversed. This will 
be important later as in this way the sign of the contri¬ 
bution to r —> i±vv can be flipped. In our notation, h is 
the SM-like Higgs, H and A are the additional CP-even 
and CP-odd Higges; due to the mainly leptophilic cou¬ 
plings of H and A, collider bounds on their masses are 
quite weak and they can be even lighter than the SM-like 
scalar boson. In particular, sin(/3 — a) = 1 always cor¬ 
responds to the SM-like limit, even for m R < rrih (this 
differs from standard 2HDM literature). 

In the following, we will assume e d = 0 for simplicity, 
as it is stringently constrained from FCNC processes [54], 
In addition, we take e u to be of the form 

/0 0 0 \ 

e“ = I 0 0 0 , (26) 

\0 x xj 

where x denotes a non-zero entry, since again e “ 3 (e 23 ) is 
severely constrained from b —> d(s )7 [91]. In the lepton 
sector we take the same structure as for the e“: 

/0 0 0 \ 

e ( = [ 0 0 0 . (27) 

\0 x xj 

In this way we avoid lepton flavour violation involving 
electrons (i.e. bounds from /1 —> ej, fi —> eee and [i —> e 
conversion in nuclei), but still allow for effects in Aa M , 
t —> hisis, and even h —► fir. 

We will not attempt to find a symmetry realisation of 
the fJ structures, but take them merely as a convenient 
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C cb > n cb 
L /°SM 


FIG. 1: Allowed regions in the C 23 -C 23 plane from 71(D) 
(blue) and 1Z(D*) (yellow) for real values of C 23 R - The lighter 
regions correspond to 2a experimental uncertainties while the 
darker regions are correspond to la. 


minimal set of parameters to explain existing anomalies. 
Should the above structures prove successful, one might 
try to find appropriate flavour symmetries to generate 
them dynamically. 


III. OBSERVABLES 


In this section we discuss the relevant processes and 
summarize the formulae needed for the phenomenological 
analysis. 


A. Tauonic B decays 


In a 2HDM of type II the charged Higgs contribution 
to B —► tv [92], 71(D) and 1Z(D*) interferes necessarily 
destructively with the SM and, in addition, 7 Z(D) and 
7Z(D*) cannot be explained simultaneously [37]. How¬ 
ever, an enhancement is possible in type-III models, see 
Ref. [51]. We do not discuss B —► tv here, which depends 
on the value of V u b- If the inclusive value is taken for V u b, 
it agrees well with the SM prediction, but it is above the 
SM value if the exclusive determination is used. (Note 
that the differences between the inclusive and exclusive 
determination cannot be explained by NP [93].) 

The SM and NP contributions relevant for 71(D) and 
7 Z(D*) are contained in the effective Hamiltonian 

'll _ zV/fr riqb . /^qb^qb . /^qb^qb 

rteff — C'SM'ASM + E R U L > 



A, H, h. H+ 


FIG. 2: Dominant one-loop contributions to r -> Ivv in the 
lepton-specific 2HDM (adapted from Ref. [65]). 


with 

OSM = T^^PLVr , 

Ot = cP R b fP L v T , (29) 

O c l = cP L b tP l v t , 


assuming massless neutrinos. The SM Wilson coefficient 
is given by Cg^ = 4 Gp V cb /i/2. The NP Wilson coeffi¬ 
cients (7fg and C£ b (given at the B meson scale), which 
parametrize the effect of NP, affect the two ratios in the 
following way [43, 94, 95]: 


K(D) 

71(D)sm 

7 Z(D*) 

7Z(D*) sm 


= 1 + 1.53? 


r C cb 


C? 1 


= 1 + 0.123? 


r^cb 

U SM 


‘ (jcb (J°b " 

R L 


1.0 


c'f 


C L 


rich 

°SM 


(30) 


r^cb 

°SM 


0.05 


(jcb _ neb 

R L 


neb 

°SM 


Here, efficiency corrections to 71(D) due to the BABAR 
detector [37] are important in the case of large contri¬ 
butions from the scalar Wilson coefficients (i.e. if 

one wants to explain 1Z(D) with destructive interference 
with the SM contribution). As shown in Ref. [50], these 
corrections can be effectively taken into account by mul¬ 
tiplying the quadratic term in C c £’[ of Eq. (30) by an 
approximate factor of 1.5 (not included in Eq. (30)). 

In this model-independent treatment, one can see from 
Fig. 1 that C|j 6 alone cannot explain 71(D) and 1Z(D*) 
simultaneously, whereas C£ b is capable of achieving this, 
e.g. with real C£ fe ~ — 1.2 |C'sm I • -^ n 0llr m °del, neglecting 
flavour-changing couplings in the down and in the lepton 
sector, only the coefficient C£ b is generated in the large 
tan ft limit 


C 


cb 



(31) 


For our phenomenological analysis we will neglect 622 and 
add the experimental (statistical and systematic) errors 
in quadrature, but include the theoretical uncertainty by 
adding it linearly on top of this. 


B. Tau decays r — > Ivv 

At tree level, only the charged scalar H + contributes 
to r -> Ivv in the lepton-specific 2HDM. (Note that the 
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100 6g 


with the loop function 


F(x ) 


1 

2 


1 + x 
4(1 — x) 


logo;. 


(37) 


In particular, F(x) < 0 (equality for x = 1), and so 
Sg < 0. In the 2HDM-X, 2 is positive and Sg negative, 
making it hard to reconcile both A e and unless rather 
high values of z ~ 0.5 are chosen, which are then in 
disagreement with the Michel parameter g [36]. 

In our model, we can however flip the sign of the r 
couplings for e 33 > m T /v, which allows for a negative z 
(Sg remains negative) as we will see in the phenomenol¬ 
ogy section. Small values z — —0.1 are then sufficient 
to satisfy the A^ constraints as well as g. Using HFAG 
values for the r decay opens up parameter space at the 
2u level, see Fig. 7. 


FIG. 3: Allowed regions for z and 5g from A M (yellow), A e 
(blue), and the Michel parameter g (red); see text for defini¬ 
tions. Darker regions are at the 2cr level, lighter regions at 
3cr one. Here we used for the experimental input the PDG 
values. 


contributions in the type-X model are the same as in the 
type-II model as the lepton sector is identical.) Due to 
the small electron Yukawa couplings, the contributions to 
Ag ree are highly suppressed compared to A[, ree , resulting 
in lepton flavour non-universality [65, 67, 68 ]. However, 
one-loop corrections can be important as well (Fig. 2), 
and contribute (to a very good approximation) univer¬ 
sally to A( [65, 67] by modifying the Wtv couplings [36] 

9wtv 9w TV (1 + 6g ). (32) 


C. Magnetic moment a M and r —> gg 


The radiative decay r —>• /ry is closely related to the 
anomalous magnetic moment of the muon. Both observ¬ 
ables are induced by penguin diagrams with internal neu¬ 
tral or charged Higgs bosons. The results can be encoded 
in the effective Hamiltonian 


"Hoff = c 


tjli 

R 


O 


R 


+ C 


ifU 

L 


o 


tfti 
L > 


(38) 


where c^' and c^ 1 are the Wilson coefficients of the 
magnetic dipole operators 


°r[ 1) = melfa^PntL/iF^ . (39) 


For <C m T and Sg <C 1, we find (notation of Eq. ( 8 )) 


A e ~ 2 Sg , (33) 

9 

TFl 

A m ~ 26g + -- - 2z—F , (34) 

4 m T 


With these conventions, the branching ratio for the ra¬ 
diative lepton decays £ t — > £/y reads 


BR(U -> l n ) 


m] 


4-7T Ti 


t(" 




J-fU 12 


)• 


(40) 


with the charged-scalar coupling (assumed to be real) 


„ _ v2 yLRH + yLRH + * 

~ ~ m2 1 W 

m H+ 


(35) 


Here we ignored flavour-changing interactions, which 
would not interfere with the SM and are tightly con¬ 
strained from flavour-changing neutral current processes. 
In addition, the H + contribution to r -> (iw leads to a 
non-zero Michel parameter g = —2z/{4: + z 2 ) [36], mea¬ 
sured to be 77 = 0.013 ± 0.020 [3]. The constraints (using 
the full expressions for Af including lepton masses [36]) 
are shown in Fig. 3 using PDG values. The SM value is 
recovered for Sg = z = 0. 

In the SM-like limit sin(/3 — a) = 1, and ignoring again 
flavour-violating couplings, we have [36] 


&9 


tan 2 /? m T 
327r 2 v 


c 33 


2 


F 


( m A 
\™ 2 h+ 


+ F 



( 36 ) 


and the contribution to the anomalous magnetic moment 
of the muon is given by 


5a y. 




(41) 


The one-loop (and dominant two-loop) Wilson coeffi¬ 
cients cl,r are given in App. A. 

It is well known that two-loop Barr-Zee-type diagrams 
[96] can dominate in some regions of parameter space due 
to enhancement factors m^/m 2 from fermions / in the 
loop, which can overcome the additional loop suppression 
ol em/tt- For the 2HDM-X, there is an additional tan (3 
enhancement for A, H , and enhanced two-loop contri¬ 
butions come from the r in the loop (Fig. 4). However, 
there are other important diagrams which could give rele¬ 
vant contributions [97 101]. While all diagrams involving 
Z or W bosons coupling to the external lepton line are 
highly suppressed, the diagrams with a closed charged 
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FIG. 4: Dominating Barr-Zee diagram of a light scalar H j? = 
A,H to (g - 2) m . 



FIG. 5: Branching ratios H + — > W + H (blue) and H + -> tv 
(red) for rriH = 30 GeV. Solid (dashed) lines are for tan /3 = 
80 (10). 


Higgs loop or a W loop connected to the external lepton 
line by a photon could be numerically relevant (replac¬ 
ing the r in Fig. 4 by a H + or a W). Without explicit 
calculation (the detailed results are given in App. A) one 
can already deduce the scaling of the relevant diagrams: 

r: rjl LR ^ tan 2 /3 (for H° = A, H), 

H+: T*° LR \ H+H - H t<mp (for H° k = H), 


E. Leptonic Higgs decays 

The decay h —> rr has been observed by CMS [102] 
and ATLAS [103] with /z-parameters (relative strength 
compared to the SM prediction) 0.78±0.27 and 1.43+(y37, 
respectively. A naive combination (also averaging the 
ATLAS result to 1.43 ± 0.40) gives 

/z TT = 0.98 ± 0.22 . (43) 


W: r R ° LR cos (a - /3) tan /3 (for H° k =H). 

Here we included only the A and H contributions, as the 
coupling of h to leptons is suppressed. The couplings of 
A to WW and H + H~ vanish due to CP conservation. 

It was pointed out in Ref. [101] that the Higgs self- 
coupling A h+h-h contribution can be very important. 
For this result, Ref. [101] allowed A h+h-h to vary be¬ 
tween —5 and 5. However, A h+h-h is not a fundamental 
parameter. It depends on the Higgs self-coupling in the 
scalar potential, but their contribution is suppressed by 
l/tan/3 and therefore negligible (see Eq. (A8)). The dia¬ 
gram involving a W loop can be important for moderate 
values of tan/3 [97], but vanishes in the SM-like limit 
sin(/3 — a) = 1. We will see later, sin(/3 — a) ^ 1 is re¬ 
quired for explaining h — > fir. Nonetheless, working at 
large tan/3, as preferred by r decays, we checked that 
also the contribution of the W diagram is numerically 
small and does not change our conclusion for (g — 2)^. 


In our model, the decay rate for h —» rr relative to the 
SM prediction takes the form 


BR ( h —> tt ) 
“ BR(h -)• rr) gM 


sin 2 a tan 2 (3 



c 33 


m T /v 


2 

(44) 


in the large tan/3 limit, assuming a SM-like production 
rate (in the SM-like limit sin(/3 — a) —> 1 we get back 
H TT —> 1). The entry e^ 2 allows for the flavour-violating 
decay h —>■ gr if cos(a — /3) ^ 0: 

BR^/zt)-^^! 2 , (45) 

where T/j ~ 4.1 MeV is the decay width in the SM for the 
125 GeV Brout-Englert-Higgs boson. 


F. Other Higgs decays 


D. Flavour changing top decays t —>• He 


Since we allow for a non-zero e^ 2 coupling, the top 
quark can decay into H k + c if the scalar H k = h,H, A 
is sufficiently light. For the branching ratio we find 


BR(f -> ff£c) 


m t 

32nTt 


r H ° LR 2 
T tc 


1 - 


m 


2 

Hi 


m 


2 

t 


(42) 


With a light H, new decay modes open up, such as 
h —> HH, A — > ZH , and H + W + H [104], followed by 
H —► rr in the large tan/3 limit. (In complete analogy 
to the case where A is light, see e.g. [105].) The decay 
rates are given by 


F (H + W + H) 


OEM m %+ ( m W \ 

48 rriyy V m H+ ) 

oem rn\ f _ m |\ 3 
48s^c^ ra| \ m 2 A ) 


(46) 

(47) 


with the top decay width T t = (2 ± 0.5) GeV. 


r(A ->• ZH) 



















7 


in the SM-like limit sin(/3 — a) ~ 1 and for rriH -C mw- 
{h -A 7777 depends on additional parameters and will 
not be discussed here.) If tan/3 is not too large, these 
channels can contribute significantly to the total decay 
rates of 77 + and A (see Fig. 5), which weakens the di¬ 
rect search bounds on these particles (these bounds often 
assume 100% decays into tau leptons for the 2HDM-X). 


IV. PHENOMENOLOGICAL ANALYSIS 

Using the formulae collected in the last section, we 
now study the phenomenology of our 2HDM and show 
that it can indeed resolve the anomalies outlined in the 
introduction. 


G. Z boson decays 


In the large tan /3 limit one expects a sizable modifica¬ 
tion of Z —>• tt at the loop level. Following Ref. [36] we 
define the ratio of decays 


_ r (z -a tt) 

T/e ~ r(Z -A ee) 


(48) 


with the experimental value 7?®* p = 1.0019 ± 0.0032 [3]. 
The deviation from the SM due to 2HDM vertex correc¬ 
tions is given by 


AR T / e 


— -^r/e 


-R 


SM 

r/e 


g|mz/67r 

T(Z —A it ) SM 


[v T 3£Au[. oop + « r SAa]. oop ] , 

(49) 


with v T = \T 3 -SwQt = Syp —1/4 and a T = \T 3 = -1/4 
being the tree level Ztt couplings and 

Av ioo P _ (^A _ 4 ^ [v T (Fi(m H ) + Fi(m A )) 

-2a T F 1 (m H +) + (v T + a T )F 2 (m H +, m H +)}, 

Ao/-°° p ^ (V “ ^ 3 ) 

-2a T F 1 (m H +) + ( v T + a T )F 2 (m H +, m H +) 

-4 a T F 2 (m H ,m A )} , 

(50) 


in the SM-like limit sin(/3 — a) = 1. The loop functions 
Fj can be found in Ref. [36]. For the region of interest 
the limits from Z decays are of similar order as those 
from r —A tvv. 

For e| 2 7^ 0, one also induces Z -A /.tr via a similar 
one-loop diagram [106]. For the mass ranges of interest 
in this work the branching ratio is approximately 

BR( Z ^).5x10-=(J444)\ (51) 

The best upper limit of 1.2 x 10 ~ 5 at 95% C.L. [3] still 
comes from LEP, but LHC searches should be able to 
improve this by a factor of few with current data [107] 
and even more with the upcoming run. 


A. TZ(D) and 7 Z(D*) 


Let us first consider 7 Z(D) and 1Z(D*). From the left 
plot of Fig. 6 we see that Ti(D) and 7 Z(D*) can be ex¬ 
plained simultaneously for negative values of e^ 2 with 
small or vanishing imaginary part. The right plot in 

„ , , , , u tan/3 (100GeV ) 2 

1 lg. 6 shows the dependence of e^ 2 on 


100 


°H+ 


requiring that F(D) and 1Z(D*) are explained. Note that 
sizable values of e^ 2 are required, i.e. of the order of 10 “ 1 . 
This is important for t —> He to be considered later. 


B. r —> ivv 

The tree-level charged Higgs contribution interferes 
destructively with the SM for eg 3 = 0. However, for 
e | 3 > m T /v the interference is constructive, allowing for 
an explanation of the PDG data, which is in more than 
a 2(j tension with the SM. The 1-loop contributions in¬ 
terfere again destructively (independently of € 33 ) and are 
important for non-degenerate A and H masses. Nonethe¬ 
less, even if rriH = 30 GeV and m A = 200 GeV, the values 
m H + = 200 GeV, e 33 = 2 m T /v and tan/3 > 60 are con¬ 
sistent with data (see Fig. 7). Furthermore, as one can 
also see from Fig. 7, for e 32 « 0.1 also 7^(77) and 1Z(D*) 
can be brought into agreement with the measurements. 
Therefore, the possibility to flip the sign of the 77, A cou¬ 
pling to taus allows us to have smaller values tan /3 than 
in the 2HDM-X. 


C. Anomalous magnetic moment of the muon 


In the anomalous magnetic moment of the muon, the 
one-loop and the two-loop Barr-Zee contribution have 
opposite sign for €33 = 0 (neglecting flavour violating 
couplings). However, for e 33 > m T /v the interference 
is constructive, allowing for an explanation with smaller 
values of tan/3 and/or heavier Higgses. Note that for 
e 33 > m T /v the 77 contribution has the same sign as the 
SM contribution while the A one has opposite sign, so 
in our scenario it is a light 77 that can solve the A 
anomaly, as opposed to a light A in the standard 2HDM- 
X. We show explicitly in the left plot of Fig. 8 that mu 
must be smaller than m A for e 33 > m T /v. As seen above, 
e 33 > m T /v is preferred by r —> [ivv. 



0.8 



FIG. 6: Left: Allowed regions in the e^-plane from 1Z(D) (blue) and 1Z(D*) (yellow) for tan/3 = 50, m H + = 200 GeV, and 
e 33 = 0. The scaling of the allowed region for e 32 with tan/3 and m H + is the same as for e 31 . e u is given at the matching scale 


m H +. Right: Allowed regions in the e 32 _ 


tan /3 (100 GeV) 


100 


l H+ 


plane for e| 3 = 0 (yellow), (blue), ^ mr (red) and ~~ ?Ut (green). 

2r^ 2r v 


D. h —> fj,r and r —> /ty 

Until now, we worked in the large tan /3 limit with 
a = 0. However, the decay h —> fj,r can only appear 
for non-zero values of cos(a — ft). In this case additional 
Barr-Zee diagrams with gauge bosons or top quarks can 
contribute to r —> /xy (see App. A). Therefore, the analy¬ 
sis is more involved than the one for the anomalous mag¬ 
netic moment of the muon. However, as we have shown 
in the case of (g — 2 ) M (where the contributions are di¬ 
rectly related to r -> /xy), cos(a — ft) 0 has actually 
only a small effect on the result. 

To explain the CMS excess in h —> /it (Eq. (10)), one 
needs a coupling strength of approximately 

sin a tan/31£ 32 1 — 3.7 x 10 -3 . (52) 

Non-zero values of e| 2 then give rise to r —> /xy. The 
experimental upper limit for r —> /x 7 is given by [108, 109] 

BR (r -7 /xy) < 4.4 x 1CT 8 , (53) 

at 90% C.L.. It is interesting to see if one can explain h —> 
/ir and (g — 2 ) M simultaneously without violating bounds 
from r —> /xy. As the loop contributions to r —» /xy 
are governed by the same Wilson coefficients as (g — 2) M 
(see Sec. IIIB), this turns out to be challenging. In the 
right plot of Fig. 8 we show the allowed regions in the 
TOa-TOjj plane for r —> /xy, h — > /xt and (g — 2 )^. As 
one can see, there is no overlap among all regions. There 
is an TOt/to^ enhanced contribution to (g — 2 )^ in the 
case e 23 7 ^ 0 and e 32 7 ^ 0. Even though we restricted 
ourselves in Eq. (27) to vanishing values of e 23 7 ^ 0, we 
checked that also for e| 2 = e 23 the effect in ( g — 2 )^ is 


small, taking into account the upper limit on e| 2 = e 23 
from r -7 /xy while aiming at an explanation of h —> fir. 

In principle, one might increase the coupling strength 
T^ A ’ H+ with the help of e 22 . This would soften the tight 
relationship 

Cr /<% * r$/T$ * 4 v/m, (54) 

originating from the dominant Barr-Zee diagram with a 
r loop which causes the incompatibility of a^ and r —> /xy 
(Fig. 4). However, a large shift in T^ A,H+ from e 22 -C 
m^/v would mean fine tuning and also strongly affect 
r -7 fxi/R. Therefore, we conclude that explaining 
and h -7 /xr simultaneously is not impossible, but rather 
difficult and would involve fine tuning. 

E. t -7 He 

For light values of m#, as preferred by the anomalous 
magnetic moment of the muon, and non-zero values of 
e 32 as required by an explanation of 1Z(D M), the flavour 
changing top decay t -7 He can have sizable branching 
ratios. In fact, as shown in Fig. 9 the branching ratio can 
be easily of the order of 10 -2 . 

V. CONCLUSIONS 

In this article we addressed the measured anoma¬ 
lies in 7 Z(D^) (3.8a), (~ 3a), r -7- /ivv (2.4a), 

and h —> /xr (2.4a) within a simple two-Higgs-doublet 
model. The Yukawa structure of our model is close 
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FIG. 7: Allowed regions in the tan /3-u /m T e^ 3 plane from 
1l(D ( * ) ) and r —t pvv at the 2 a level. The yellow region is 
allowed by r —> rvv using the HFAG result for ran = 30 GeV 
and rriA = 200 GeV, while the (darker) blue one is the allowed 
region using the PDG result. The red, orange, green and 
magenta bands correspond to the allowed regions by 
for different values of £ 32 - The gray region is excluded by 
Z — > tt and r — > evv. For mn ~ rriA the allowed regions 
from r —» ri/v would be slightly larger. 


mH + =200Gev 



to the lepton-specific 2HDM (type X), but with some 
additional Yukawa couplings involving third-generation 
fermions that give rise to the b-c (necessary for 1Z(D^)) 
and fi-r transitions (relevant for h —¥ rt) as well as cor¬ 
rections to tt couplings (important for r —► rvv). 

Let us summarize the logic of the article. 

• r —» rvv prefers e | 3 > m T /v. If one wants to 
account for the PDG average, also tan/3 > 50 is 
desirable. 

• ci/j, favours small values of m# for e 33 > m T /v as in 
this case the Barr-Zee contribution with a r loop 
has the correct sign and the diagram involving a 
Higgs self-coupling can be relevant. 

• 7 Z(D) and 1Z(D*) point towards quite large nega¬ 
tive values of e 32 (of order — 0 . 1 ). 

• In case of a light H (as preferred by a^), sizable 
decay rates for t —>• He are possible if 7 Z(D) and 
7 Z(D*) are explained via e 32 . This decay could be 
observed at the LHC in the process pp —)• it, t —>■ 
He, H — > TT. 

• h —► /ir can be explained using e| 2 . In this case a 7 ^ 
0 and constraints from r —> /zy arise. As the Barr- 
Zee contributions in r —> /zy are directly correlated 
to the ones in a^, a simultaneous explanation is 
difficult. 


• If one attempts to explain h — > rt (disregarding 
a^), the exotic process pp —► it, t — 1 He, H —► /it 
can occur at the LHC. 

Therefore, the future prospects are very promising: the 
decay h —> /it implies rates for r — > /zy in reach of fu¬ 
ture experiments. More data on tau leptons is neces¬ 
sary to test our model, in particular r —► tvv, r —► /zy, 
and h —I tt. Furthermore, the light H and the flavour- 
changing couplings required for IZ(D^) lead to the de¬ 
cay t —► He, followed by H —> tt (or even H — > /it). 
which can be searched for at the LHC. While we did not 
attempt to find a symmetry realisation of the pattern as¬ 
sumed for the e{j structures, it would be very interesting 
to find appropriate flavour symmetries to generate them 
dynamically, as the model works very well phenomeno¬ 
logically. An additional venue of interest would be the 
inclusion of dark matter in order to explain the galactic 
centre gamma-ray excess [ 110 ]. 
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Appendix A: Wilson coefficients for a M and li —> if 7 


The effective Hamiltonian relevant for ii —► f/y and 
is given in Eq. (38) with operators from Eq. (39). The 
Hermiticity of the Hamiltonian allows us to deduce cr 
from cr in complete generality via 


ifii _ m l f Zdf* 

Cr - Cp 

me. H 


(Al) 


so we will only show cr in the following. The final Cl.r 
of course requires a sum over all the individual Cl } r we 
present here. 

At one loop the neutral Higgs (F7° = H, h, A) penguin 
contribution to * is given by 


tfU 

C RH° 


3 

V _—_ 

. . , 1927r 2 m^ 0 
k,j =1 M k 


r H° LR H° k LR * 
1 ifij 1 titj 


J_' 

I L f P r 

me- 3 1 


me j 

me i 


r H° LR r H° LR 

1 £ftj 1 ijii 


^9 + 6 In 



(A2) 
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FIG. 8: Left: Allowed regions in the mA-rriH plane from the anomalous magnetic moment of the muon at the 2 a level for 
tan/3 = 70, rn H + = 200GeV, cos(a — /?) = 0 and different values of £ 33 . Right: Allowed regions in the mA-rriH plane from 
(g — 2) m , r —>■ g-y and h — > gr for tan(/3) = 30 and £33 = 2 m T /v. For h —> r^i blue corresponds to cos(cr — (5) = 0.1 and light 
blue to cos(q — /3) = 0.2. The allowed region for Aa M is maximal in the sense that we allowed for the three possibilities £32 7 ^ 0, 
£32 = £23 7 ^ 0 and £32 = —£23 7 ^ 0 as the latter ones can give m T /m ^ enhanced one-loop contributions. However, the effects 
turns out to be small, as £ 32,32 is stringently constrained from r —> p 7 . In addition, we checked that the effect of A h is very 
small. 


and the charged Higgs penguin contribution takes the form 



FIG. 9: The contour lines denote BR(t —> He) x 100 as a 
function of £32 and rriH■ The coloured regions are allowed by 
7 Z(D) and 1Z(D*) for different values of tan /3/m 2 H+ . 


3 

t-ft-i _ e r ff+ LR r H+ LR * 

RH+ 3847r 2 m 2 „ + m t . ^ ^ 

H+ j=l 


(A3) 


It is well known that two-loop Barr-Zee-type diagrams 
[96] (see Fig. 10) can dominate in some regions of pa¬ 
rameter space due to enhancement factors m 2 /m 2 from 
fermions / in the loop, which can overcome the additional 
loop suppression ctem/tt- For the 2HDM-X, there is an 
additional tan /3 enhancement for A, H , and enhanced 2- 
loop contributions come from the r in the loop. However, 
there are in addition other important diagrams [97-101] 
that could give relevant contributions to be summarized 
and converted to our conventions in the following. 


1. Diagram (a) 


The most relevant Barr-Zee diagram contains an ad¬ 
ditional charged fermion / in the loop (Fig. 10a). 


e f £i,(a) _ 


128tt 4 


N fQ 2 f 


mi. 


E 

i=h,H,A 


■pz pi 


mj 9i ^ 


(A4) 


where Nj and Qf are the number of colours and charge 
of fermion /, respectively, and r l j = m 2 /m 2 . The loop 
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FIG. 10: Barr-Zee diagrams relevant for the Wilson coefficients c 


l,r ■ 


function is given by 
9i a \i")= [ dx 

Jo 


W.W (A5) 


jo x i l -x)-r 

where 

N h (x) = N h (x) = 2x(l - x) - 1, Na(x) = -1 

2. Diagram (b) 


(A6) 


The diagram from Fig. 10b contains a charged scalar 
in the loop and hence depends on the H + H~H 3 self¬ 
interactions. The relevant dimensionless trilinear Higgs 
couplings are given by 


1 


A H+H-h — -2 — m h) 


2 \ sin(2a + 2/3) sin(a — /3) 


\ mH+ 2 


sin(4/3) 
sin(a — /3) 


+ AA tan(2/3) cos(o: + /3), 

(A7) 

1 / 2 2 \ sin(2(t T 2/3) cos(cr /3) 

A h+h-h = ym H -m h ) - 


sin(4/3) 


1 


m. 


„ | m~„+ -) cos(a — /3) 

+ AA tan(2/3) sin(a + (5 ), 


(A8) 


and \ H+H ~ A = 0. AA = A 2 —Ai is a free parameter [111]. 
While the last two terms in Eq. (A8) vanish in the SM- 
like limit sin(/3 — a) = 1, the first is still suppressed by 
l/tan/3 for large tan/3. The Wilson coefficient is then 
given by [101] 


C R 


128 v / 2tt 4 m ti 

-pi 

„ ST" 1 e f e * 


i—h.H,A 


C l ^H+H-H a 9i 


( b) 


m 


H+ 


(A9) 


where = —( H = —( A = 1 and the loop function 


is 


&z( r ) = / d;c 


x(l — x) , / r 


In 


x(l — x) — r \x(l — x) 


(A10) 


3. Diagram (c) 

The contribution of the diagram in Fig. 10c with a W 
in the loop is [101]: 


JfUM _ G F e 3 


E 


647T 4 . t—* me. 

i—h.H 


-r 


(All) 
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where 1Zhi = sin(/3 — a), IZhi = — cos(/3 — a) and the 
loop function is 


9h%( r ) = \f o da ' 


x(l — x) — xr[3x(4x — 1) + 10] 


5(1 — x) — 


x In 


; (1 — x ) 


(A12) 


4. Diagrams (d), (e), (f) 


Finally, there are three types of Barr-Zee diagrams, 
where the virtual Hf and 7 propagators are replaced by 


H + and W: Figs. lOd (where the fermion line is a t/b- 
or 6/t-loop), lOe and lOf. Defining the loop function 

f r a x+r b (l-x) \ 

G(r a , rb ) = -77 - ^ t~ X) 1 -T . ( A13 ) 

x(l — x) — r a x — r°(l — x) 

and the matrix 

sin(/3 — a) cos(/? — a) 0\ 

TZ = | — cos(/3 — a) sin(/3 — a) 0 j , (A14) 

0 0 1 / 


the diagrams in Fig. lOd, lOe and lOf are [101] 
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(A17) 
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